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ABSTRACT 

We performed numerical simulations of star cluster encounters with Hernquist's 
TREECODE in a CRAY YMP-2E computer. We used different initial conditions (rel- 
ative positions and velocities, cluster sizes, masses and concentration degrees) with 
total number of particles per simulation ranging from 4608 to 20480. Long term inter- 
action stages (up to 1 Gyr) when the pair coalesces into a single cluster are compared 
with isolated LMC clusters. Evidence is found that when seen in a favourable plane 
these resulting clusters show elliptical shapes due to the disruption of one of the com- 
panions. These elliptical shapes are essentially time independent but they do depend 
on the initial structural parameters of the pair components. We also analyzed the 
fraction of stars that are ejected to the field by the interaction. Wc found that this 
fraction can be almost 50% for the disrupted cluster. These simulations can represent 
a possible mechanism to explain the ellipticity observed on several star clusters in the 
Magellanic Clouds. 
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1 INTRODUCTION 

The study of star cluster pairs in the Magellanic Clouds 
has been explored extensively in the past years (Bhatia 
& Hatzidimitriou 1988; Bhatia & McGillivray 1988; Bica, 
Claria & Dottori 1992; Rodrigues et al. 1994, hereafter 
RRSDB; Bica & Schmitt 1995; de Oliveira et al. 1998, here- 
after Paper I; Bica et al. 1999). The evolution of physical 
pairs can provide fundamental insight into the past history 
of cluster formation in the Magellanic Clouds. 

In the last decade, N-body simulations of stellar sys- 
tem encounters have been the main tool to investigate the 
dynamical processes that can occur in such interactions, 
like mergers and tidal disruption. Many of the first simula- 
tions worked with equal mass encounters (White 1978; Lin 
& Tremaine 1983; Barnes 1988). Different mass encounters 
have been carried out by Rao, Ramamani & AUadin (1987), 
Barnes & Hut (1986, 1989), RRSDB and Paper I. These 
studies indicate that tidal disruption and merger are two 
important processes in the dynamical evolution of a binary 
stellar system. However, there remain important issues not 
yet well tested, as the long term stages of such interactions, 
and the fraction and how stars are stripped from the cluster 
by the parent galaxy tidal field. 

Some isolated clusters in the Magellanic Clouds present 
interesting structures with higher ellipticities than Galactic 
globular clusters (van den Bergh & Morbey 1984). Several 
explanations have been proposed like the presence of sub 
clumps merging with the main cluster which could produce 



the impression of high ellipticities in these clusters (Elson 
1991). Another possibility is that old mergers in advanced 
evolutionary stages could show such ellipticities (Sugimoto 
& Makino 1989). 

This work is a continuation of Paper I where we ana- 
lyzed the morphology of selected cluster pairs in the Mag- 
ellanic Clouds and have compared them to those obtained 
from numerical simulations of star cluster encounters. A pre- 
liminary discussion was given in de Oliveira et al. (1999). In 
the present paper we study the long term behaviour of some 
of our numerical models (up to 1 Gyr), and compare the fi- 
nal stages of these simulations (isodensity maps, ellipticities, 
isophotal twisting) with the structure of isolated clusters. 

In Section 2 we describe the method and the conditions 
employed in the present simulations. In Section 3 we describe 
the Magellanic Cloud cluster images and the procedure to 
derive the isodensity maps and ellipticity measures. The nu- 
merical results and discussions are presented in Section 4. 
In Section 5 we give the main conclusions. 



2 THE METHOD AND THE INITIAL 
CONDITIONS 

We performed the simulations using TREECODE (Herquist 
1987) in the CRAY YMP-2E computer of Centre Nacional 
de Supercomputacjao of the Universidade Federal do Rio 
Grande do Sul (CESUP-UFRGS). For a complete descrip- 
tion of the method see Paper I. 
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Table 1. Input parameters for the simulations, by columns: (1) 
number of particles of the cluster; (2) cluster total mass; (3) half- 
mass radius; (4) maximum radius of the cluster; (5) mean velocity 
of the particles (modulus) and (5) softening parameter. 



Table 2. Collision parameters. By columns: (1) elliptic orbit 
model; (2) total number of particles of the two clusters; (3) mass 
ratio; (4) eccentricity of the orbit; (5) pcricentre of the orbit; (6) 
initial relative velocity . 



Model 


Npart 


Mt 


Rh 


J^max 


ymd 


t 






(Mo) 


(pc) 


(pc) 


(Km/s) 


(pc) 


A 


16384 


10^ 


4.94 


20.0 


3.53 


0.30 


B 


4096 


10* 


3.72 


15.0 


1.94 


0.47 


C 


512 


10^ 


2.00 


8.0 


0.83 


0.50 


D 


512 


103 


4.00 


8.0 


0.62 


0.50 



Model 


Npart 


hh 


/hh 


ro{pc) 


e 


p{pc) 


V, (Km/s) 


E9AB10 


20480 




10 


70.0 


0.9 


20.0 


4 


E9BC10 


4608 




10 


37.0 


0.9 


10.0 


2 


E9BD10 


4608 




10 


37.0 


0.9 


10.0 


1 


E7BD10 


4608 




10 


37.0 


0.7 


10.0 


1 


E6BD10 


4608 




10 


37.0 


0.6 


10.0 


1 



2.1 The Initial Conditions 

In Paper I, cluster N-body models were generated for 16384, 
4096 and 512 equal mass stars, corresponding to total masses 
of 10'' M©, 10'' Mq and 10^ Mq, respectively. The cutoflF radii 
of the Plummer clusters are 20, 15 and 8 parsecs, respec- 
tively. These values are comparable to diameters found in 
Magellanic Clouds clusters (Bhatia et al. 1991). We gener- 
ated two different concentration degrees for the 512 particle 
cluster, one with half-mass radius 1/2 of the cutoff radius 
and another with 1/4 (Table hi). The interaction models are 
characterized by a pericentre distance p (distance of closest 
approach) and the initial relative velocity Vi. These initial 
velocities are comparable to the observed random velocities 
in the LMC disk (Freeman, Illingworth & Oemler 1983). The 
initial conditions (t=0) used in Paper I were: (i) the large 
cluster (in the case of equal mass clusters, the more con- 
centrated one) is located at coordinates (X,Y) = (0,0); (ii) 
the small cluster lies at a distance ro = 5 x rt (where rt 
is the tidal radius). Beyond this distance the tidal effects 
are negligible. The initial relative velocity at this distance 
was obtained from the two-body formulae. The position and 
velocities of the particles during the encounters were com- 
puted in the centre-of-mass frame of the entire system. At 
various times, the essential data containing the positions 
and velocities of all the particles were stored for later anal- 
ysis. The computation was stopped when disruption of the 
cluster occurred or, in the case of open orbits, when the rel- 
ative separation of the clusters reached distance d > 2ro. 
We point out that the softening parameter adopted for the 
interaction is the smallest value of the two cluster models 
(column 7 of Table 1). Indeed, with reference to the AB clus- 
ter interaction (which has the largest e difference) we carried 
out a simulation where cluster B is allowed to evolve isolated 
adopting the cluster A softening value for At ~200 Myr (~ 
2 relaxion times). Analysis of the time evolution of cluster 
B structure shows that it is essentially the same for e = 0.30 
as that obtained for the cluster with e = 0.47. Consequently 
in the AB cluster interaction significant structural changes 
are not expected by the adoption of the smaller e value. 

In the present Paper, we proceeded the simulation from 
this point on for some models with closed orbit where dis- 
ruption occurred (see Table B) . We used the same procedures 
as in Paper I. These simulations were ran up to 990 Myr for 
the models presented here. 

We illustrate in Fig. hi the time evolution of an ellip- 
tic orbit encounter involving two clusters with 16384 and 
4096 particles respectively (model E9AB10, see Table 0). 



In this Fig. we can see the small cluster being deformed by 
the massive one with the formation of a bridge (~50 Myr), 
subsequently complete disruption occurs (~135 Myr); at the 
end of this simulation the clusters coalesce into a single one, 
with the stars of the disrupted cluster forming a halo around 
the final cluster, and some of them being ejected. 



3 ISODENSITIES FOR LMC ELLIPTICAL 
CLUSTERS 

In the final stages of our simulations the pairs coalesced 
into a single cluster with a distinct structure as compared 
to the original ones. We noticed that in viewing angles close 
to the original orbital plane of encounter, the final single 
cluster presented an ellipticity larger than those of the ini- 
tial clusters (our models have a spherical symmetry). In this 
paper we have selected some LMC clusters with morpholo- 
gies resembling those of the present models (Sect. 4.4). We 
checked many clusters with important ellipticity, as indi- 
cated in previous studies (Geisler & Hodge 1980, Zepka & 
Dottori 1987, Kontizas et al. 1989, 1990). In particular the 
selected clusters present increasing ellipticity outwards. In 
Table H we show these clusters with their V magnitudes, 
SWB types derived from interpreted UBV photometry (Bica 
et al. 1996, Girardi & Bertelli 1998) and corresponding 
ages. Also there are available ages from colour-magnitude 
diagrams: (i) NGC1783 has 700- 1100 Myr depending on 
adopted distance modulus (Mould et al. 1989). Based on 
the same data, using the AV turnoff/clump method, Geisler 
et al. (1997) obtained 1300 Myr; (ii) NGC1831 has an age 
of 400 Myr according to Hodge (1984). CCD data provided 
500-700 Myr (Mateo 1988) and 350-550 Myr (Vallenari et 
al. 1992), depending on adopted models; (iii) NGC2156 has 
age ~60 Myr according to Hodge (1983); (iv) NGC1978 has 
an age ~2000 Myr (Olszwewski 1984, Bomans et al. 1995). 
Geisler et al. (1997) derived 2000 Myr by means of the /W 
turnoff/clump method. There is good agreement of ages de- 
rived from the CMD studies and the integrated colours (Ta- 
ble ||). NGC1783 and NGC1831 have ages comparable to the 
long term stages of our models, while NGC2156 falls some- 
what short, and NGC1978 has an age ~ a factor 2 larger. 

The images of this selection were obtained from the Dig- 
itized Sky Survey (DSS). The plates are from the SERC 
Southern Sky Survey and include Ilia- J long (3600s), V 
band medium (1200s) and V band short (300s) exposures. 
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Figure 1. Time evolution of an encounter (model E9BD10) projected in the XY plane. 



Table 3. Ages for the clusters in the present sample. By columns: 
(1) cluster name; (2) observed visual magnitude; (3) SWB type; 
(4) age. 



Name 


V 


SWB 


Age 






type 


(Myr) 


NGC1783 


10.93 


V 


900 


NGC1831 


11.18 


IVA 


400 


NGC1978 


10.70 


VI 


3000 


NGC2156 


11.38 


III 


140 



The PDS pixel values correspond to photographic density 
measures from the original plates, and are not calibrated. 

The digitized images, similary to those generated by 
the models, were treated with the IRAFn package at the 
Institute de Fi'sica — UFRGS, applying a 2-d Gaussian filter 
to smoothen out individual stars, and to create isodensity 
maps. The ellipticity measurement was made with the IRAF 



* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation, U.S.A. 
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package task Ellipse which fits elliptical isophotes to data 
images with the use of a Fourier expansion; 



J/ = ^ A„sin{nE) + B„cos{nE) 



(1) 



where E is the position angle. The amplitudes 
{A3, B3,A4,B4,), divided by the semi-major axis length and 
local intensity gradient, measure the isophotal deviations 
from perfect ellipse. Note that the TRAP routine requires 
that one indicates a first guess for the object major axis 
position angle. Subsequently the routine iteractively fits the 
best solution for each isophote. We stored the coefficient Bi 
which tells whether the ellipse is disk (a positive B4, param- 
eter) or box-shaped (a negative B4 parameter), Bender et 
al. 1989. We also stored the semi-major axis position angle 
variation. It should be noted that in general ellipticity does 
not present simple behaviour. Ellipticity varies with radius 
(Zepka & Dottori 1987, Kontizas et el. 1989,1990) and so it 
is difficult to assign one to a cluster. We measured ellipti- 
cal isophotes starting at Rh and stopping at the last possible 
non diverging ellipse. In the models, this last ellipse occurred 
around the radius containing 90% of the cluster mass. 



4 DISCUSSION 

We considered models in which the less massive cluster (per- 
turbed) is allowed to move in elliptic orbits around the more 
massive one (perturber). The perturbed cluster is assumed 
to move in an orbit of eccentricity e = 0.6, 0.7, 0.9. The 
values of the pericentre and the eccentricity are meaningful 
only if the clusters are assumed to be point masses moving 
in Keplerian orbits. However, soft potential orbits are not 
conical and our use of these definitions is meaningful only 
to a good approximation, but not strictly (White 1978). 

The collision parameters of the simulations are given 
in Table H. The model designation (column 1 of Table Cl) 
contains information on the encounter conditions. E refers 
to an elliptic orbit encounter; the number following the letter 
E is related to the orbital eccentricity (column 5 of Table m); 
A, B, C and D refer to the model type (Table hi); the last 
number is the orbit pericentre (column 6 of Table H). 

From the analysis of the time evolution of each model 
as illustrated in Fig. hi for E9AB10, it is possible to observe 
the trend to form bridges in the beginning of the encoun- 
ters (see also Paper I). As time goes by the smaller cluster 
is completely disrupted by the massive one, resulting in a 
single cluster at the end of the simulation. 



4.1 The Structure of the Final Stages 

In Pig. H, we illustrate the same encounter as in Fig. hi but 
in a different plane projection (parallel to the orbital plane). 
We can observe the disruption of the smaller cluster, with 
its stars occupying orbits around the massive cluster. As 
these stars remain preferentially orbiting in the same plane 
of the original encounter, they introduce an elliptical shape 
to the final single cluster, when seen in a favourable plane 
projection. 

In Pig. H we show isopleths, ellipticity, position angle 
and 4*'' cosine coefficient (hereafter B4 ) for the model 



E9AB10 in a XZ projection. These measures were made 
for four different evolution times, the first panel being the 
initial massive cluster A before interacting (t=0 Myr). We 
can observe a radial variation of ellipticity during the inter- 
action, which has also been observed for all other models. It 
also can be seen a general trend of decreasing ellipticity to- 
ward the inner parts of the models. Zepka & Dottori (1987) 
and Kontizas et al. (1989) have observed internal variations 
of ellipticity in a number of LMC clusters, finding a gen- 
eral tendency of increasing ellipticity towards the inner parts 
of these clusters. However, they reported a few exceptions 
(~ 5% of the data) which in turn have counterpart in our 
models, and could be the result of an interaction. Indeed 
there are about 300 cluster pairs among a total population 
of 7847 extended objects in the SMC, inter-Cloud region and 
LMC (Bica & Schmitt 1995 and Bica et al. 1998). At least 
50% of the pairs are expected to be interacting (Bhatia & 
Hatzidimitriou 1988) and we need a favourable projection 
(nearly edge-on orbits) to observe elliptical structure like 
those of the models. 

In Fig. y there is little time variation of ellipticity during 
the lifetime of our model. This variation occurs mainly in 
the outer parts of the cluster, probably by the continuous 
escaping stars of the disrupted cluster. 

There occur small position angle variations (Fig. 0), 
which indicate that the elliptical shape of the cluster stands 
out in the plane of the interacting encounter. The B4, coeffi- 
cient clearly has a positive peak, which is a sign of the pres- 
ence of a disc component structure in the resulting model. 
This suggests that the elliptical shape in our final clusters 
is due to a rotation disk formed mainly by stars of the dis- 
rupted cluster. 

In order to check whether mass or size could affect the 
behaviour of our models, we did the same analyses for mod- 
els with different mass and concentration degree, but with 
same orbital encounter. In Fig. we present the same mea- 
sures of Fig. for other two models at the final stage and 
we see similar ellipticity variations in radius. In order to 
compare absolute ellipticity variations between models, and 
also between different time stages in the same model, it is 
required to define ellipticity at a common radial distance. 
As we have the main shape changes due to variations in 
the outermost parts of these clusters we decided to measure 
the ellipticity as the mean between 3Rh and the outermost 
ellipse measured. 

In Table we give the mean ellipticity between 3Rh and 
the outermost ellipse Bout for all the models in different time 
stages. We conclude that in model E9AB10 Cout shows a 
small tendency of increasing with time; differences among 
the models are found mainly in the shape of e versus radial 
distance curve. 

When comparing ellipticity between different models at 
the same appropriate time (Table m, we observe a small re- 
duction of Cout from a more massive model (E9AB10) to 
a less massive model (E9BC10). This relation (ellipticity) 
versus (mass) has been suggested for the LMC clusters (van 
den Bergh & Morbey 1984) where high mass clusters have 
higher angular momentum or have more difficulty in shed- 
ding it than do low mass clusters. However, when comparing 
a more massive model (E9AB10) with a less massive one 
(E9BD10) — where we have in the last a different con- 
centration degree for the disrupted cluster — we observe no 
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Figure 2. Time evolution of an encounter (model E9BD10) projected in the XZ plane. 



significant variation in Cout- Thus, the initial concentration 
degree affects the final shape of the resulting cluster. This 
suggests that a more concentraded cluster has a deeper po- 
tential well, making more difficult stripping its stars by the 
more massive cluster. On the other hand, stripping stars in 
the less concentrated cluster is more efficient, resulting in 
a more pronounced halo expansion. This is supported by 
the results in Table 5. This halo expansion when seen in an 
edge-on view seems to be an important factor to stablish the 
shape of the final merger. 

We also compared the possible ellipticity variation be- 
tween models with equal mass and concentration degree, but 



different initial orbital parameters. We observed no signifi- 
cant differences in Bout between the models. 



4.2 Rotational Velocity Field 

In Fig n we show the line-of-sight rotational velocity field 
(Vy) for the model E9AB10 along the major axis at 990 Myr 
seen in a XZ projection. The merger consists approximately 
of a rigid rotation core and an outer halo ( r > 20 pc) with a 
Keplerian fall. The rotational velocity along the major axis 
has a peak value of about 1 Kms~^ at a radius of 20 — 
30pc. It is interesting to see that the rigid rotation extends 
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Figure 3. In each panel we have the isopleth, the eUipticity, position angle and Bi coef. for the model E9AB10, in a XZ projection. 
These measures were made for four different evolution time as indicated. The first stage shows the initial massive cluster A before 
interacting. 

Table 4. Mean eUipticity measured between 3 Rh and the outermost ellipse for the models in different time scales. By columns: (1) 
model; (2) to (9) eUipticity for different evolutionary times and respective mean errors. 



Model 


time 

(Myr) 


e 


terror 


time 

(Myr) 


e 


terror 


time 

(Myr) 


e 


€-error 


E9AB10 


225 


0.12 


0.008 


500 


0.13 


0.005 


990 


0.15 


0.008 


E9BC10 


- 


- 


- 


500 


0.10 


0.01 


990 


0.08 


0.009 


E9BD10 


- 


- 


- 


500 


0.14 


0.01 


990 


0.14 


0.004 


E7BD10 


- 


- 


- 


500 


0.14 


0.01 


990 


0.13 


0.009 


E6BD10 


- 


- 


- 


500 


0.16 


0.006 


990 


0.14 


0.004 
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Figure 4. Same as Fig. M for otiier two models at the final stage (990 Myr). 



exactly to the radius of the initial massive cluster before 
the encounter. This indicates that the main contribution to 
the velocity field for r > 20 pc is due to stars of the small 
cluster. This becomes evident when we plot the two star 
cluster member sets for a merged final model (Fig. pi) in 
a XZ projection. In this Fig. we plot the model EABIO at 
T=990 Myr and a reference box (side of 40 pc) centred at the 
more massive cluster. It is clear that the main contribution 
to the final composite system for r > 20 pc is due to stars of 
the small cluster. Velocity distributions obtained with high 
resolution spectra with large telescopes might reveal such 
rotation curves, while in turn could be a signature of evolved 
stages of cluster merging. 



4.3 Mass Loss 

In Table H we give the radii containing 10%, 50% and 90% of 
the mass for the massive cluster in each interaction (model 
A or B) for f = and for the final cluster stage. We can 
observe that Ri,o varies little for the initial and final clus- 
ter stage when comparing models with same initial mass, 
suggesting that the massive one has its original stellar con- 
tent little affected by the collision. However model E9AB10 
shows a contraction of the core, showing that the initial mass 
may affects the cluster final structure. The radius containing 
90% of the mass is larger for all our models, which shows 
that the final cluster outer halo is formed mainly by stars 
from the disrupted cluster. The latter result is in agreement 
with those presented in Paper I for the radial mass distribu- 
tion. There we showed that the maximum expansion for the 
disrupted cluster in E models always occurs in the plane of 
the encounter and is a minimum in the z direction. When 
we compare the projected density distribution between the 



Figure 5. The linc-of-sight rotational velocity field along the 
major axis of the model E9AB10 at 990 Myr. The cluster is seen 
in a "edge-on" projection (XZ plane) where the abscissa is the X 
axis. 



models before interacting and at the end of the simulation 
we see that they are similar within the region r < 10 pc 
(Fig m). However, the final cluster extends with a depen- 
dence of p ~ r~'^. 

Together with changes in the internal structure of our 
final single cluster, the present simulations suggests that a 
fraction of cluster stars may be ejected to the field due to 
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Figure 6. For the merger EABIO at T=990 Myr we separate the 
stars belonging to each cluster. The left panel shows the stars 
of the more massive cluster, and the right one stars of the less 
massive in a XZ plane projection. We also plot a 40 pc side box 
centred at the more massive cluster as a reference. 



Figure 7. Projected density distributions for the model E9AB10 
at 990 Myr (dashed curve) and for the model A before the inter- 
action (solid curve) . The radius is in parsecs and the density is in 

MQ.pc~^. 



Table 5. Radii containing 10%, 50% and 90% of the mass before 
and at the final stage of the encounters (in parsecs) . 



Table 6. Final cluster escaping mass fraction. By columns; (1) 
model; (2) mass fraction (3) absolute mass loss. 



Model 




RiO 


^50 


R90 


A (t=0) 




2.03 


4.94 


11.79 


B (t=0) 




1.55 


3.71 


8.89 


E9AB10 (t= 


=990 Myr) 


1.00 


3.44 


17.62 


E9BC10 (t= 


=990 Myr) 


1.34 


3.89 


12.83 


E9BD10 (t= 


=990 Myr) 


1.33 


3.79 


13.77 


E7BD10 (t= 


=990 Myr) 


1.35 


3.86 


13.94 


E6BD10 (t= 


=990 Myr) 


1.29 


4.01 


14.13 



the encounter. In order to quantify this fraction, we calcu- 
lated the total energy per particle, identifying those that 
had positive total energy which are expected to escape the 
system. 

In Table H we have for each model the final fraction 
of stars that escape the cluster (in units of solar mass). The 
mass fraction that is ejected to the field is very small, consid- 
ering the ensemble of the pair (< 3% ), so the contribution 
to the field stars by a cluster pair encounter is not very sig- 
nificant. It should be noted that when we study the merger 
of the cluster pair we are not taking into account the tidal 
field of the parent galaxy. This should truncate the outer 
halo of the final cluster, which would contribute with an 
increasing fraction of stars ejected to the field. Adopting a 
tidal radius rtidai ~ 65pc (distance moduli of 18.5 for LMC, 
Westerlund 1990) for observed cluster's halo (Elson et al, 
1987), we can estimate how much mass is beyond this limit. 
For the model E9AB10, we estimate M ^^ 50% of the mass 
of the disrupted cluster (4.5% of the total mass) is beyond 
this limit. So, it can be concluded that although the total 



Model 


mass fraction 


absolute mass 




(%) 


(Mq) 


E9AB10 


2.90 


3195 


E9BC10 


2.88 


316 


E9BD10 


3.20 


352 


E7BD10 


1.85 


204 


E6BD10 


1.15 


127 



mass loss for the pair is not very significant, the assumption 
of a maximum radius for the final cluster, implies that half 
of the disrupted cluster stars feed the field. 



4.4 Comparison of Simulations with Magellanic 
Cloud Clusters 

Isopleth maps of projected planes at a given time f of a 
suitable model can be compared to the observed isodensity 
maps of selected Magellanic Cloud clusters to infer on their 
dynamics. In Paper I we have used this method to search for 
evidences of interacting pairs. In the present study, we use 
model isopleths together with other measured parameters in 
order to compare our simulation final stages with possible 
observational evolved products of cluster pair interaction. 
Some examples are given in Figs. 0, H, FUj and [111 where 
we show some isolated LMC clusters isophotes, and radial 
dependence of ellipticity, position angle and B4, coefficient. 
The examples shown here have an ellipticity radial vari- 
ation with a trend to increase towards the outer parts. NGC 
1783 has an ellipticity curve compatible with our model 
E9BD10 when seen in a XZ projection at 990 Myr (Fig. 0, 
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right). The cluster age (Sect. 3) is fully compatible with the 
model evolutionary time. This scenario can be explained by 
the result of a cluster pair encounter with the subsequent 
merger of the pair in the early history of NGC1783. 

NGC1831 has many resemblances with the models. The 
ellipticity (Fig. 0) increases outwards, but in the exter- 
nal isophote the value decreases again. This is observed in 
the model E9BC10 at 990 Myr (Fig. |, left). The cluster 
presents a considerable isophotal twisting like in the model. 

As pointed out in Sect. 3, the age of NGC2156 falls short 
of the long term stages of the models. But notice that asym- 
metrical early stages (Fig. g) at ~ 25 Myr , if seen favourably 
in a given direction might present ellipticity variations. The 
negative B4 coefficient behaviour suggests that if the cluster 
is indeed a merger, disruption has not yet occurred to create 
a disc shape. 

NGC1978 is an interesting intermediate age cluster due 
to the pronounced ellipticity (Fig. hlj), which is already 
presented in the innermost available isophote. No isopho- 
tal twisting occurs in this case. Recently Kravtsov (1999) 
searched for special variations in the colour-magnitude dia- 
gram of NGC 1978. They found evidence of some variations 
which they attributed to a possible metallicity spread. Bo- 
mans et al. (1995) did not find any evidence of age variation 
within the cluster. A merger scenario does not require dif- 
ferences in age and/or metallicity. Indeed two clusters born 
in the same association are expected to have a close orbit 
encounter, which is the more favourable case for interactions 
(Paper I). 

The Bi coefficient has difi'erent behaviours in the 4 clus- 
ters, and in some cases has varying values in a given clus- 
ter. We might be witnessing clusters with either discoidal 
(positive values) or boxy (negative values) shapes. The av- 
erage value of B4, for NGC1783 suggests an overall boxy 
shape. NGC2156 is also boxy, but with a systematic trend 
to become disc-shaped in the outer parts. NGC1831 shows 
a predominant positive B4, suggesting a disc-shape. Finally, 
NGC1978 on the average is positive thus indicating a pre- 
dominantly disc-shape. However the cluster has two distinct 
B4 value regions, 20< r <50 is definitively disc while r > 50 
should be classified as boxy. Behaviours like this have coun- 
terparts in elliptical galaxies (see e.g. the sample of Goud- 
frooij et al. 1994). The B4, coefficient is a promising tool to 
explore the possibility of merger in galaxies or star clusters. 



Figure 8. In each panel we have the isophote, the ellipticity, the 
position angle and the 64 coef. for the cluster NGC 1783. 



5 CONCLUSIONS 

We used N-body simulations to study the final morphology 
and structure of star cluster pair encounters. We also com- 
pared these morphologies with those of elliptical isolated 
LMC clusters. The main conclusions may be summarized as 
follows: 

(i) Close orbit encounters of cluster pairs can lead to a 
single final cluster with a distinct structure of the original 
ones. When seen in a favourable plane projection, models 
show an elliptical shape comparable to those of some ob- 
served isolated LMC clusters. This suggests that tidal en- 
counters could be a mechanism to explain the ellipticity of 
several clusters in the Magellanic Clouds. 

(ii) Evolved stages appear to be stable after ~ 200 Myr, 
suggesting that the resulting ellipticity is not transient. The 



Figure 9. In each panel we have the isophote, the ellipticity, the 
position angle and the S4 coef. for the cluster NGC 1831. 



simulations indicate that the initial concentration degree 
and mass affect the final shape of the resulting cluster. In 
the final stages all models present halo expansion, and some 
present core contraction. 

(iii) The models show mass loss for the composite sys- 
tem.The fraction of stars ejected to the field by the en- 
counter, is not so significant (~ 3% with respect to the sum 
of both clusters) . However it can represent up to ~ 50% of 
the stars of the disrupted cluster, if we assume a tidal radius 
of Ttidai ~ 65 pc for the final merged cluster. 
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Figure 10. In each panel we have the isophote, the ellipticity, 
the position angle and the B4 coef. for the cluster NGC 2156. 



Figure 11. In each panel we have the isophote, the ellipticity, 
the position angle and the B4 coef. for the cluster NGC 1978. 



(iv) The velocity distribution of some final stage models 
presents characteristic velocity patterns in favourable plane 
projections, which could be used as an observational con- 
straint to test the present scenarios. 

Finally, we call attention that the merging of spheroidal 
systems can produce disc-shaped products. Thus not only 
boxy-shaped systems might be related to mergers. 
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